Binary mixture of pseudo-spin-^ Bose gases with interspecies spin exchange: from 
classical fixed points and ground states to quantum ground states 
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We consider the effective spin Hamiltonian describing a mixture of two species of pseudo-spin-^ 
Bose gases with interspecies spin exchange. First we analyze the stability of the fixed points of 
the corresponding classical dynamics, of which the signature is found in quantum dynamics with 
a disentangled initial state. Focusing on the case without an external potential, we find all the 
ground states by taking into account quantum fluctuations around the classical ground state in 
each parameter regime. The nature of entanglement and its relation with classical bifurcation is 
investigated. When the total spins of the two species are unequal, the maximal entanglement at 
the parameter point of classical bifurcation is possessed by the excited state corresponding to the 
classical fixed point which bifurcates, rather than by the ground state. 
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I. INTRODUCTION 

A remarkable discovery in recent years is that bifurca- 
tion in classical dynamics is related to quantum entangle- 
ment in the ground state of the corresponding quantum 
Hamiltonian. In addition to its theoretical demonstra- 
tion in the Dickc model [3, HI, in a model of two cou- 
pled giant spins describing magnetic clusters p|, and in 
an integrable dimer model Q, this association has also 
been studied in the of Bosc-Einstein condensation 
(BEC), including two-component BEC [H, @ and two- 
mode atomic- molecular BEC @, HI ■ More recently, a 
classical bifurcation has been observed in an experiment 
realizing an internal Josephson effect in a spinor Bosc- 
Einstcin condensate, as an important step toward entan- 
glement generation close to critical points 0] . Moreover, 
in a laser-cooled atom, experimental evidence has been 
observed for entanglement being a quantum signature of 
chaos [l(J. 

On the other hand, a novel kind of BEC, the so-called 
EBEC, that is, BEC of interspecies spin singlet pairs, 
was found to be the ground state of a Bose system com- 
posed of two species of pseudo-spin- ^ Bose atoms with 
both intraspecics and interspecies spin-exchange interac- 
tions in a considerable parameter regime |llT[l5| . Under 
the usual single orbital-mode approximation, the Hamil- 
tonian of this system can be transformed into that of 
two coupled giant spins. Alas, the ground states in all 
parameter regimes have not yet worked out. 

In this paper, we make each of the above two lines 
of research useful for the other. We focus on the case 
in the absence of an external potential. First, we study 
the bifurcation of the classical dynamics corresponding 
to the Hamiltonian of this Bose mixture, by analyzing 
the stability of each fixed point. Quantum dynamics 
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displays some features similar to the classical dynamics 
if the initial state is a disentangled state, which, how- 
ever, is not an energy eigenstate. When the numbers of 
the atoms of the two species are equal, a bifurcation of 
the fixed points indeed corresponds to a quantum phase 
transition to a maximally entangled ground state. When 
they are unequal, the quantum state corresponding to 
the classical fixed point which bifurcates is also maxi- 
mally entangled at the bifurcation point. However, it 
is not the ground state. Finally, we analytically obtain 
all the quantum ground states by considering quantum 
fluctuations around the classical ground state in each pa- 
rameter regime. The analytical results fit the numerical 
results very well. 

The rest of this paper is organized as follows. The 
model is introduced in Sec. |TIJ The fixed points and bi- 
furcations are studied in Sec. IIII1 with the detailed cal- 
culation reported in the Appendix. The classical ground 
states are described in Sec. lIVl and the classical evolution 
is described in Sec. El whose quantum analog is described 
in Sec. I VII Section IVIII describcs the absence of the corre- 
spondence between classical bifurcation and maximal en- 
tanglement in the case of unequal populations of the two 
species. In Sec. IVIII1 we find the quantum ground state 
in each parameter regime by approximating the Hamil- 
tonian around the classical ground state there. Finally 
the paper is summarized in Sec. IIXI 



II. THE MODEL 

Consider a dilute gas composed of two distinct species 
of Bose atoms with the following property [l2|, • Each 
atom has an internal degree of freedom represented as 
a pseudospin- i , with z-component basis states f and J.. 
Under the usual single orbital-mode approximation, for 
each species a(a = a, b) and pseudospin a (a =t, 40; oru y 
the single-particle orbital ground state </w(r) is occu- 
pied, then the Hamiltonian can be transformed into that 
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of two coupled giant spins with spin quantum numbers 
S a = Ng/2 and Sb = Here we focus on the uniform 

case [1J, [la , for which 



From the Hamiltonian ([I]), one obtains the equations 
of motion 



7~L — J±(SaxSb x + SaySby) + J z S az Sb z 



(1) 



where J± = 4irh 2 t; e / (m a b^) while J z = 47r?i 2 (£ s — 
£d)/( TO ab^), with m a i, being the reduced mass of an a 
atom and a b atom, £ e is the scattering length for the 
scattering in which an a atom and a b atom exchange 
pscudospins, £ s being the scattering length for the for- 
ward scattering in which an a atom and a b atom have 
different pscudospins without spin exchange, £j being the 
scattering length for the forward scattering in which an 
a atom and a b atom have the same pseudospin with- 
out spin exchange, Q being the volume of the system, 
S a = aJ-Sovciv , s aa i is the single spin operator, a a de- 
notes the annihilation operator associated with <p aa (r) of 
species a. Without loss of generality, suppose S a > Sb- 

The corresponding classical Hamiltonian is obtained 
from ([!} by treating the spin operators as the classical 
spin variables, 

7~L c i = J±(S ax Sbx + S a ySby) + J z S az Sbz (2) 
= Jxy/ffi - Sl z ){Sl - Si) cos(^ a - <p b ) 
+J z S az Sb z (3) 



dS a 



dS a 



dt 
dt 



J J_S aZ S fjy J Z S a ySp Zl 

= ~ J ±S az S p x + J Z S aX S/3 Z , 

= J±(S a ySf3 X S aX Sf)y) . 



(4) 



The corresponding classical equations of motion, ob- 
tained either from the classical Hamiltonian ([2]) or from 
the quantum equations of motion by replacing the spin 
operators as the classical spin variables, can be written 



as 



dA 
~dt 



JA, 



(5) 



where A — (S ax , S a y, S az , Sb x , Sby, Sb z )^ 
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In studying the stability of a fixed point, J becomes the 
Jacobian matrix when the spin variables adopt the values 
at this fixed point. 

The classical Hamiltonian in the form of (|3]) suggests 
that the classical state of the system is completely deter- 
mined by the variables S az , Sb z , and ip a — <~Pb- 

We shall use © studying the evolution in Sec.fV] while 
using ([2]) in analyzing the fixed points in Scc. lIIIl because 
there is arbitrariness of angles <j> a and 4>b in some fixed 
points. 

III. FIXED POINTS AND BIFURCATIONS IN 
CLASSICAL DYNAMICS 

The fixed points of the classical dynamics are obtained 
from 



The stability of each fixed point can be examined first 
by studying the eigenvalues of Jacobian matrix J at this 
point: It is stable if every eigenvalue has a negative real 
part, while it is unstable if any eigenvalue has a positive 
real part. Otherwise, one cannot judge whether the fixed 
point is stable from the eigenvalues of J, but it is cer- 
tainly stable if a Lyapunov function can be constructed. 
A Lyapunov function T is such that in a neighborhood 
of the fixed point, C is minimal (or maximal) at the fixed 
point, and dC/dt < 0(or dC/dt > 0). 

There exist eight fixed points in our system. Their 
stability is analyzed in the Appendix. The stable pa- 
rameter regimes of these fixed points are shown in 
FigQ] We specify the fixed point in terms of direc- 
tion n a of the spin vector S Q = S a n a , that is, n Q = 
(sin 8 a cos ip a , sin 9 a sin ip a , cos 9 a ) , with < 9 a < 7r,0 < 
(p a < 2tt, a — a,b. 

As shown in FIG. [TJ the fixed points and their stable 
regimes are the following. 
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FIG. 1: (Color online) The stable parameter regimes of 
the fixed points on J± — J z plane. The fixed points are 
(I) n a = -n b = (0,0,1) or n a = n b = (0,0,-1), (2) 
n a = -n b = (cos tp, sin tp,0), (3) n a = n b = (cos tp, sin tp, 0), 
(4) n a = n b = (0,0,1) or n a = n b = (0,0,-1), (5) n a = n 6 , 
(6) n a = — n b , (7) n a = (sin#cos(/3,sin#sin</5,cos#), 
while rij, = (sin #cos tp, sin 6 sin tp, — cos 8), (8) 
n a = (sin 8 cos tp, sin 8 sin tp, cos 8), while n;, = 
(— sin 8 cos tp, — sin 6 sin tp, cos 0). The fixed point (4) is 
always stable, which the stable regimes of the other fixed 
points are indicated by using the curves with arrows at both 
ends. 



(1) n a = nb = (0,0, ±1). One spin is the parallel 
to the z direction while the other is antiparallel to the z 
direction. This fixed point is stable when 771 1 J z \ > \J±\, 

where 771 = \ (^/f + ^/f) . 

(2) n Q = — rib = (cos tp, sin tp, 0), where < tp < 2ir. 
The two spins are antiparallel and arc both on the x — y 
plane. This fixed point is stable if Jl > and Jl > 772 J 2 , 
or Jj_ < and Jj_ < 772 Jz, where 772 = i'^ . 

(3) n a = rift = (cos tp, simp, 0). The two spins are par- 
allel and are on the x — y plane. This fixed point is 
stable when J± > and J± > —rfrJz, or J± < and 
Jl < -mJz- 

(4) n a = rij, = (0,0, ±1). The two spins are both 
parallel or antiparallel to the z direction. This fixed point 
is always stable. 

(5) n a = — n b . The two spins are always antiparallel. 
This fixed point only exists at Jl = J z and is stable. 

(6) n a = nb. The two spins are always parallel. This 
fixed point only exists at Jj_ = J z and is stable. 

(7) n a = (sin0cos</5,sin0siniy9,cos0) while rib = 
(sin 6 cos tp, sin 6 sin tp, — cos 9). The z components of the 
two spins are opposite. This fixed point only exists at 
Jl = — J z and is stable. 

(8) n a = (sin^cost/JjSin^sin^jCOsfl) while nb = 
(— sm6 cos tp, — sin # sin (p, cos 6). The xy components of 
the two spins are opposite. This fixed point only exists 
at J_l = — Jr and is stable. 




FIG. 2: (Color online) Classical ground states in different 
regimes of parameters J± and J z . They are just the eight fixed 
points, but there are differences in the parameter regimes al- 
though there are overlaps, (i) J z > \J±\, the fixed point 
(1) ; (ii) J± > I J* I, the fixed point (2); (iii) J± < -\J Z \, 
the fixed point (3); (iv) J z < —\J±\, the fixed point (4); 
(v)J± = J z > 0, the fixed point (5); (vi) J± — J z < 0, 
the fixed point (6); (vii) J± — —J z > 0, the fixed point (7); 
(viii) J± — ~J Z < 0, the fixed point (8). 



It can be seen that J± = ±771 J z and J± = ±772 J z are 
bifurcation points. 



IV. CLASSICAL GROUND STATES 

As depicted in Fig. [31 it can be found that classically 
the energy is minimal at fixed point (1) when J z > | Jj_|; 
at fixed point (2) when J± > \J Z \; at fixed point (3) when 
Ji < — |J Z |; at fixed point (4) when J z < — | Jj_|; at fixed 
point (5) when J± = J z > 0; at fixed point (6) when 
J_L = J z < 0; at fixed point (7) when J± = —J z > 0; at 
fixed point (8) when J± = — J z < 0. 

If S a = Sb, we have r\\ = 772 = 1; therefore the param- 
eter regimes of the bifurcation points are completely the 
same as those of the classical ground states, respectively. 
But if S a 7^ Sb, there are differences although there are 
overlap regimes. 



V. CLASSICAL EVOLUTION 

Because the Hamiltonian conserves S az + Sb z , we study 
the dynamical evolution of S az — Sb z and (tp a — ipb) /2 for 
some given values of S az + Sb z - When S az = S a while 
Sbz = Sb, or S az — —S a while Sbz = —Sb, the system 
is at the fixed point (1). When S az = 0, Sbz = while 
(tp a — tpb)/2 = tt/2, the system is at the fixed point (2). 
When S az = 0, Sbz = while {tp a — ipb)/2 — or n, the 
system is at the fixed point (3). 

We have studied the evolution trajectories near fixed 
points (1), (2) and (3) under various values of S a = Sb- 
As shown in Fig. [31 when J±/J z < 1, fixed points (1) 
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FIG. 3: (Color online) The evolution trajectories on the plane 
of (Saz — Sbz)/Sb and (ip a — fb)/2, for various values of S a = 
S b . Here r, = J x /J„ S az + S b z = 0. 
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FIG. 4: (Color online) The evolution trajectories on the plane 
of (Saz — Sbz)/Sb and (ip a — v'b) /2, for various values of S a = 
2S b - Here r\ = J±/J z . The solid lines describe the dynamics 
for Saz + Sbz = and the dot lines describe the dynamics for 

Saz + Sbz = ±(S a — Sb)- 



and (3) are stable while fixed point (2) is unstable; when 
J±/Jz > 1, fixed points (2) and (3) are stable while fixed 
point (1) is unstable. This conclusion is reached by con- 
sidering that a fixed point is stable if the evolution tra- 
jectories are loops around a fixed point, otherwise it is 
unstable. 

We have also studied the case of S a ^ S b . As shown 
in Fig. H] for S a = 2Sb, the evolution trajectories are 
different from the case of S a = S b - For J±/J z = 0.9 and 
for J±/J z = 1.03, three fixed points are all stable. 

Note that all the results of the numerical simulation 
are consistent with the above theoretical analysis. 



VI. QUANTUM EVOLUTION WITH 
DISENTANGLED INITIAL STATE 

To simulate a quantum process closest to classical evo- 
lution, we choose as the initial state a disentangled state, 



FIG. 5: (Color online) Quantum dynamics of {S a z)/S a for 
S a = Sb — 300 and various values of r\ = J± / J, ■ The figures 
on the upper line exhibit the stability of fixed point (1). The 
figures on the second line exhibit the stability of fixed point 
(2). The unit of t is 1/J Z . 



which can always be written as 



0) = (e-W 2 cos^| t>« + eWW£| l) a ) N -> 



)(e -^/2 cos || t)b + e W2 siri || ;)b) ^ 

= \S a n a ) ® \S b n b ). 



(7) 



where | ^) a denotes the spin state of a single atom of 
species cv, while \S a n a ) represents the state of the to- 
tal spin of species a. Tn this state, the spin compo- 
nents of each species are similar to classical spins; that 
is, (S ax ) = S a smd a cosip a , (Say) = S a sm0 a smip a , and 
(Saz) = S a cos0 a (a = a, b). Moreover, we choose a and 
(p a in such a way that (S Q ) corresponds to a fixed point 
in classical dynamics. For fixed point (1), n a = — n b = 
(0,0, 1); thus the initial state is \ip {1) ) = j t)t N "\ i)f Nb - 
For fixed point (2), n a = —rib = (cost/9, sin 0); thus 
the initial state is |V> (2 )) = (^~ ilp/2 \ t>« + 7f c ^ /2 | I 
)„)® JV »(^e- < ^ 2 | f>6 + 7se'W2| y b )»N K This is so be- 

cause Saz = N a ^ — ^f- = — N a _i, where N aC r is the 
number of atoms of species a with spin a (a = a, 6; 

<7=t4)- 

For each initial state in the form of (J7J, we evaluate 



(Saz(t)) = (i;\e lHt Saze- lHt \iP), 



(8) 



whose evolution actually represents the change of the dis- 
tribution of the atoms of species a in the two pscudospin 
states. 

We choose the same initial conditions as in the classi- 
cal case in last section to start the quantum dynamics. 
It is found that the classification of the stability of the 
classical dynamics still applies. The result is shown in 
Fig. [5] for the case of S a = Sb, and in Fig. [5] for the case 
ofS a = 2S b . 
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FIG. 6: (Color online) Quantum dynamics of {S a z)/S a for 
S a = 2Sb = 300 and various values of rj = J±/J z . The figures 
on the upper line exhibit the stability of fixed point (1). The 
figures on the second line exhibit the stability of fixed point 
(2). The unit oft is l/J z . 



FIG. 7: (Color online) The numerical result of the entangle- 
ment entropy of the quantum states of S z — and S z = 
S a — Sb, as a function of r\ . Here S a = 3Sb, J z > 0. 



The reason why quantum dynamics under the disen- 
tangled initial state is so close to classical one is the fol- 
lowing. In Hciscnbcrg picture, the quantum equations of 
motion Q reduce to the classical ones JS]), with (S a i), 
(i = x, y, z), substituting the corresponding classical spin 
variable. 



VII. BIFURCATION AND ENTANGLEMENT 

The ground state can always be written as 
\ G s z ) = fi m i S z )\S a ,m) a \S b , S z - m)b, where 
the interspecies entanglement can be quantified as 
-E m f 2 (™,S z )\og 2Sb+1 f 2 (m,S z ) 0. It has been 
shown that when S a = Sj, = S, the entanglement of the 
ground state is maximal at J± = J z , where the ground 
state is {V2STT)- 1 J2 s m= _ s (~l) m \S,m) a \S,-m)b 
[13]. Using the transformation U = e l7TSb ~ , we 
can obtain the ground state at J± = —J z as 

(y/2S + l)" 1 Yt, m =-s 1^, m) a \Sb, -m)b, which is 
also maximally entangled. 

When S a = Sb, there is only one bifurcation point at 
J± = J z between the fixed points (1) and (2); there is 
also only one bifurcation point at J± = —J z between the 
fixed points (1) and (3). Each of these bifurcation points 
corresponds to a maximally entangled quantum ground 
state. 

However, when S a ^ Sb, r\\ ^ r]2, there are two bi- 
furcation points Ji = rji J z and J± = r/2J z between the 
fixed points (1) and (2). Similarly, there are two bifur- 
cation points Jj_ = —r)iJ z and J± = —r\iJz between the 
fixed points (1) and (3). 

Numerical calculations indicate that in consistency 
with the classical ground states, the total z-component 
spin exhibits the following features. When J z > while 
7? < 1, S z = ±{S a - S b ). When J z > while T) > 1, 
S z = p, with p = if S a — Sb is an integer while p = ±1/2 



if S a — Sb is a half integer. 

Numerical results of the entanglement entropy of the 
states with S z = S a — Sb and S z = 0, varying with 77, 
are shown in Fig [7] for some integer values of Sb and 
S a = 3Sb- For 77 < 1, the ground state is the one with 
S z = S a — S b , whose entanglement values are plotted 
as empty triangles. For 77 > 1, the ground state is the 
one with S z = 0, whose entanglement values are plotted 
as filled triangles. Therefore, there is a discontinuity of 
entanglement in passing 77 = 1, where both states are 
degenerate ground states. 

According to the bifurcation analysis discussed above, 
when S a = 3S b , the two bifurcation points are 771 = 
±(y<T/3 + a/3) w 1.1547 and 772 = 0.6. As indicated 
in FIG[71 the entanglement entropy of the states with 
S z = and S z = S a — Sb is maximal at 771 and 772, re- 
spectively, and decreases rapidly in deviating from each 
of them, with the decrease more rapid for 77 larger than 
the maximal point. 

Therefore, when S a 7^ S b , the quantum state corre- 
sponding to each classical fixed point still possesses max- 
imal entanglement at the parameter point where the fixed 
point bifurcates. However, this quantum state is not the 
quantum ground state. In other words, the entanglement 
of the quantum ground state at each bifurcation point is 
not maximal anymore. 



VIII. ANALYTICAL SOLUTIONS OF THE 
QUANTUM GROUND STATES 

We now proceed to analytically find out the quan- 
tum ground states on all Jj_ — J z parameter regimes, 
by using effective Hamiltonians which describe devia- 
tions from the classical ground state in each parameter 
regime. All the ground states are summarized in Fig. HI 
Regimes A (J z > |Jj_|) and B (J z < — |Jj_|) both cor- 
respond to |£ s — £d\ > |£ e |; i.e., the interspecies spin ex- 
change scattering is quite weak, with A and B differing 
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FIG. 8: (Color online) The quantum ground states in all 
the parameter regimes. \Gi) (i — A, B, C, D) represent the 
ground states in the four bulk regimes, and are given in Sees. 
IVIH Al to I VIII Dl (1) to (4) represent the four boundaries, in 
which the ground states are described in Sec. IVIH El 



n a = 0, 1, 2, • ■ • , 2S a . When S a is very large, (f£f a ) < 

2S a , S a - PS (2S a ) 1 ^ 2 f al S az Sbz ~ Sbfafa + S a f b + fb — 

S a Sb- Then the Hamiltonian JT]) can be approximated as 



H PS - J z S a S b + USbflfa + S a f b fb) 
+ J±V^S~b(fVt + fbfa), 

Then we make the Bogoliubov transformation 



fc 



Ai + 1 



f a + sgn(J±) 



O J b ' 



f d = sgn(J±) 



o J a * 



Ai + l 



fb, 



(12) 



(13) 



where sgn(J±) is the sign of Jj_, Ai = 
—, J *( s *+ S ») . when S a = S b and J z = J ± , 

Ai ± 1 should be 1. Hamiltonian ([12]) becomes 



in whether the equal-spin forward scattering length is 
larger or smaller than the unequal-spin forward scatter- 
ing length. Regimes C (J± > \J Z \) and D (J± < —\J Z \) 
both correspond to |£ e | > \£ a — £<j|; i-e., the interspecies 
spin exchange scattering is quite strong, with C and D 
differing in whether the spin-exchange scattering length 
is positive or negative. 

Regime A (J z > J±) corresponds to £ s — £d > |£e|> i-e., 
the interspecies spin exchange scattering is quite weak. 

A. J z > \J ± \ 

As shown in Fig. 2, in this parameter regime, the clas- 
sical ground state is fixed point (1), i.e., \S a , S a ) \Sb, — Sb) 
or \S a , -S a )\S b , S b ). 

First we consider the quantum ground state near 
\S a , S a )\Sb, —Sb). One can make the Holstein-Primarkoff 
transformation |16j |. 

Sa— = f a^^S a ~ fa fa, S a + = \J 2S a — fafafa, 

S a z S a ft fa, . . 

/ / ^ ' 

Sb- = \/2Sb - flf b f b , S b+ = ft \j2Sb - ft A, 

Sbz — f b fb - Sb, 

with S a ± = S ax ^iS a y, a = a,b, f a and f a being bosonic 
operators satisfying f a \n a ) = ^/n^\n a - 1), fl\n a ) = 
y/n a + l\n a + 1), [f a ,fp] = 0, [f a , fj\ = 5 a p, where 

\n a ) = \S a ,S a ~n a ) a , (10) 
\n b ) = \S b ,-S b + n b )b, (11) 



H A = e lc ftfc + eidftfd + E 10 , (14) 

where E w = - J z S a S b + id^ll {Sa + S b ) - 
\J±\y/ (Aj-l)S7s; , eic = J ^- 1] S a + J ^ +1) S b - 
|Jx|V (A^-l)S^ e ld = J - (A 2 1+1) ^ a + J ' (A 2 1 - 1) S 6 - 
| J±\y/(Af — l)S a Sb- Thus the energy spectrum is 
E A (n c ,n d ) = ei c n c + e ld n d + E w , where n c and n d are 
nonnegative integer numbers. For J z > \J±\, eic and e\ d 
are always positive; therefore the ground-state energy is 
Si (0,0) = E w . When J z -> \J±\, E o (0,0) approaches 
— Sb(S a + 1), which is the the exact ground-state energy 
at Jz = \J±\- 

Like the original Hamilton ([TJ, the effective Hamilto- 
nian (fT2"j) also conserves the z component of the total 
spin. Therefore any of its eigenstates can be written as 

\ipi(n c ,n d )) = >J g\{n c , rig, m)\S a , m) a \S b , S z - m) b , 

ra 

(15) 

where max(— S a ,S z — Sb) < m < 

min(5' a ,5'2 "b Sb), 

gi(n c ,n d ,m) is the expansion coefficient, and S z is the 
total z component of the spin system. Using (|13|) and 
considering that \ipi(n c , n d )) is an eigenstate of both /J/ c 
and f\fd, with eigenvalues n c and n d respectively, we ob- 
tain 

n c -n d = S a - Sb- S z . (16) 

For the ground state |'0i(O, 0)), S z = S a — S b . 

It is easy to find the ground state \ipi(0,Q)) of (fl"2|) 
from/ c ^i(0,0))=0, 



7 



|^(0,0)) = D J2 " s 3"(^) 

m=S a —2Sb _ 

s b 

= D —sgn(J± 



Ai + 1 



m— — Sh 



Sa-Sb— m 



\S a , m) a \S b , S a — Sb — m)b, 



Ai + 1 
Ai-1 



|S ,S a — Sb — m) a \Sb,m)b, 



(17) 
(18) 



where D = [(|^)^ ^gS^^ ]' 1/2 is 
the normalization coefficient. When J± — > 0, |Ga) — > 
\S a , S a }\Sb,—Sb), which is an exact ground state of the 
Hamiltonian ((T|) with J± = and J z > 0. 

The excited states of (fl~2|) can be obtained by the action 



of /t and f d on the ground state 0)). With 5 a > 

Sfc, ei c < eid, for a given S^, the lowest excited state 
is |^i(n c ,0)) if S z < S a - S b and is |^i(0,n rf )) if S z > 
S a — Sb- These two excited states can be written as 



Oa-nc I a -. 

|ViK,0)> = D c eM-K*S az )J2 E (-ir^V x^r +2n C c ^_ m |5 ,m> |5 6) 5' a -5' 6 -n c -m) b , 

ii=0m=S a -2S b -n 1 

(19) 



|^(0,n d ))=£) d exp(-iC7r5 az ) £ (- J ^±i)™^/c^ a+nd _ m \S a , m) a \S b , S a -S b + n d - m) b , (20) 

m=S a -2S b +n d V 1 



where Q — if Jj_ > while £ = 1 if Jj_ < 0; D c and 
are the normalization constants, and G™ is the binomial 
coefficient. 

Now we consider the ground state close to the other 
classical degenerate ground state \S a , — S a )\Sb, Sb), in a 
way similar to the above. The Holstcin-PrimarkofF trans- 
formation is 



S' b _ = J2S b - f'Jfl S' b+ = J2S b - fff'bfi 



Sb 



fb fbi 



S a - = \J%S a — fa fafai S a ^ 
S az = fa fa ~ Sa: 



faX/^Sa f a f a , 



where the bosonic operators f' a and f' h act on 



\n a 
In! 



\S a: S a 



L b j = \Sb, Sb — n b )b- 
Thus one obtains a Hamiltonian 

H' A f» —J z S a Sb + J z (S b fa fa + Safb fb)' 



(21) 



(22) 
(23) 



eiJ'Jf' c + e ld f'Jf' d + E 1(u 



where 



f'c 



Ai + 1 . . / Ai-1 -, t 

O fa + S 9 n (Jx)\ o fb > 



f> i = a gn(J ± )J^-±% 



Ai +1 



(24) 



ft 



Therefore the eigenstates can be written as 
Wi{ n 'a n 'd)) = ^29i(n' c ,n / d ,m)\Sa,m) a \S b , S z - m) b , 

m 

(25) 

with the constraint 



n 'c - n 'd = Sa - S b + S z 



For the ground state |^i(0, 0)), S z = S b - S a . 



(26) 
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St 



K(o,o)) = ^' 



m— — Sb 



-sgn(J±) 



Ai-1 



\S a ,S b - S a - m) a \S b ,m) 



(27) 



where D' 



_ r (Ai+i) 2 



-(Ai-l) 2 



2(A 2 -l) 6 i. 



-]- 1 ' 2 . When J ± 



0, \Ga) — > \S a , — S a }\Sb, Sb), which is an exact ground 
state of the Hamiltonian (JT]) with J± — and J z > 0. 

The excited states of (f!M]) can be obtained by the action 
of f'J and f'J on the ground state \ip{(0,0)). With S a > 
Sb, tic < eidi f° r a given S z , the lowest excited state is 
KK, 0)) if 5, > S b - S a and is |^(0, n' d )) if 5 Z < S b - 
S a . The explicit expressions of \t/j[(n c , 0)) and |^((0, n<j)) 
are like dHJ) and (gUJ) for |^i(n c ,0)) and \tpi(0, n d )), with 
|5 a , m) , \S b , S a - S b - n c - m) b and |5 6 , S a -S b + n d - 
m) replaced as |5 a , — m) a , \S b , —S a + S b + n c + m) b and 
\Sb, — S a + Sb — n d + m), respectively. 

It is important to note that |-0i (0, 0)) and |'0i(O, 0)) are 
orthogonal unless S a = S b . Therefore, when S a ^ Sb, the 
ground states are doubly degenerate ones |^i(0,0)) and 
1^(0,0)) at each parameter point in this regime. 
When S a = S b , 7 = (^i(0,0)|^(0,0)) = 
2(25 i ,+i)(A 1 -i) . j^nce we mus ^ g nc j the ground 



(A 1 + l) 2s b + 1 -(A 1 -l) 2s i. + 1 

state in their two-dimensional subspacc. Clearly 
(Vi (0,0)|H|^(0,0)) = (^(0,0)1^1^(0,0)) « E w . 
(Vi(0,0)|H|^i(0,0)) = (^(0,0)|H|Vi(0,0)) w e 1o7 . 
Consequently, the ground state is found to be 



\G A (S a = S b )} 



V2 



M(0,0)) + K(0,0)>], (28) 



with energy £ 10 (1 + 7). The energy of ^[|V'i(0,0)) - 

M(0,0))] is ^0(1-7). 

When J± and J z approach the boundary J z = —J± > 
from the regime of \Ga), \Ga) approaches e l7rSaz \S a — 
Sb,±(S a — Sb))- When J± and J z approach the boundary 
J z = J± > from the regime of \Ga), \Ga) approaches 
\S a — Sb,±(S a — Sb))- 



B. J z < -\Jj_\ 

In this parameter regime, the classical ground states 
are \S a , S a )\Sb, S b ) , in which the two spins are both along 
the z direction, and |iS a , —S a )\S bl —Sb), in which the two 
spins are both along the —z direction. 



Consider the ground state close to \S a , S a ) a \S b , S b ) b . 
We make the Holstcin-Primarkoff transformation S a - = 



h-a\f 2S a hah a , S a + — "\J2S a hah a h a , S az — S a 

h' a h a , where h a and h' a are bosonic operators, now with 
\n a ) = \S a , S a — n a ) a , where n a = 0, 1 • • • , 2S a . When 
S a is very large, (h+h a ) « 2S a , S a - w {2S a ) x / 2 hl, 
S a zS bz « S a Sb—Sbh a ~h a —S a h b hb- Then the Hamiltonian 
(UJ becomes 



Hi 



~-J z S a S b - J z (S b h) a h a + S a h\h b ) 



+ J±y/S^b(hih b +hlh a ). 
We define another two bosonic operators, 



(29) 



; i-a 2 £ / 1 + A 2 , 

h c = -sgn(J±)\l — - — h a + \l — - — h b , 



h d = sgn(J±)\j 1 + „^ 2 h a 



where A 2 = ^ (5a ~ 5b) , 

ytj2(S a -S b ) 2 +4J 2 ± S a S b 

nian (UJ) becomes 



1-A 2 



(30) 



Then the Hamilto- 



Hb = ^2ch\h c + t2dh d hd + E20, 



(31) 



where E 2 



JzS a S b , £2c 



r J,(l+A 2 ) 



S a 



Ml - A2) S b + \J ± W(l-Al)S a S b l e 2d = - l h{1 - A2) S a - 



JA1 + A2) S b - \J±_W{1- A 2 2 )S a S b }. The energy spectrum 
is E B {n c ,n d ) = e 2c n c + e 2d n d + E 2 o- Thus the ground 
state is \S a , S a )\S b , S b ). 

The excited state \ip2(n c ,n d )) of (|29)l can be ob- 
tained by the action of h\ and h' d on the ground state 
\S a , S a )\S b , S b )- It is obvious that 62c is always larger 
than e 2 d, for a given S z , the lowest excited state is 
|^2(ra c ,0)), with 



Sa ~t~ S b S z 



(32) 



h - a / — 

\Mn c ,0)) = D 2c eM~^S az ) (~ \j Y^^y" ^ C nJS a , S a ~ m) a \S b , S b - n c + m) b} (33) 

m— max(n c — 2Sb ,0) 

i 

where C = if J ± > while C = 1 if J± < 0, and D 2c = [ (i+A2)"-(i-^-a+^ 
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where no — max(n c — 2Sb,0). 

Using the same method, we obtain the approximate 
Hamiltonian H' B close to the other classical ground state 
\S a ,—S a )\Sb,—Sb), which turns out to be the quantum 
ground state. The Holstein-Primarkoff transformation is 



now S a+ = h'£y2S a - h'ah'a, S a - = \j2S a - h'£h' a h' a , 

S az = h'^h' a — S a , where h' a and h'£ are bosonic opera- 
tors, with \n' a ) = \S a , n' a 

— S a )cf Thus we have Eqs. (|29[) 
to (f3lj) with primed operators. In this set of eigenstates, 
the lowest excited states for a given S z are \ijj 2 {n' c ,0)), 
with 



n — S z + S a + Si, 



(34) 



\ip 2 (n c , 0)) are like ([33]) . with 15a, S a — m) a and \Sb, Sb — 
n c + m)b replaced by \S a , ~S a + m) a and \S b , —Sb + n c — 
m) b . 

1^2(0,0)) and 1-02 (0, 0)) are orthogonal, hence are just 
the doubly degenerate ground states in this regime, and 
can be written as 



\G B ) = \S a + Sb,±(S a + S b )). 



(35) 



When Jl and J z approach the boundary J z = J± < 
from the regime of \Gb), \Gb) approaches \S a + 
Sb,±(S a + Sb)}- When J± and J z approach the boundary 
Jj_ = — J z > from the regime of \Gb), \Gb) approaches 
\S a + S b ,±(S a + Sb)). 



C. Jx > I J* I 

In this parameter regime, it is convenient to rewrite 
the Hamiltonian as 



H = Jxy/(S* - S% Z )(SI - Si) cos(^ a - Vb ) + J z S az S bz , 

(36) 

where (p a (a = a, b) is the azimuthal angle. 

In the vicinity of the classical ground state, S az ~ 0, 
Sb z ~ , (fi a — fb ~ I", for simplicity we define ip a — ipb = 
Wab + 7T. Therefore 

-He « -J±y/S a S b {Sa + l)(Sb + 1) + ~ i ± ~ J \ ; Si 



J±£,+ - Jz 



2z 



4( - Jz) 

J±.t- a V2 



J±Z+ - J:. 



■S u 



+ 2J ± ^S a S b (S a + l)(S b + l)(^ff, 



(37) 



where S 2z = S az -S b z, = H^+f a ), £- = Uk'fJ- 
S z commutes with H, and is thus a constant of motion. 
Then P 2 = S 2z ~ 7^177^ and X 2 = <^ a b/2 arc conju- 
gate variables, as S az and tp a are canonically conjugate 
variables. The Hamiltonian is then similar to that of a 



harmonic oscillator. The energy spectrum is thus 



E 3 (n,S z ) = -J ± y/S a S b (S a + l)(S b + l) 

72 _ 72 

, J ± J z C 2 



4( J±Z+ - Jz) z 

+ (n + i)^/2J ± (J ± e+ - Jz)y/S a (S a + l)S b (S b + 1) 

M - Jj_ ^^(^ + 1)^ + 1) + .,'l\~ Jz ^ S't 



4(J±£+ - Jj.) 



+ (n + i)V2J±(J±C+ - Waft, 



(38) 



where n is the quantum number of the harmonic oscilla- 
tor. The eigenstate for n = can be written as 

1^3(0, S z )) = Z(S Z ) ^ /( TO > S z )\S a ,m) a \S b , S z - m)b, , 



(39) 



where f(m,S z ) = (— l)™ 1 cxp[— 



2J± SaSb 



min(S a ,S z + S b ); Z(S Z 

ization coefficient. 

The ground state is thus 



S z )], with max(-S a ,S z - Sb) < m < 

is the normal- 



\Gc) = |^3 (0,p)) = Z(p)^2f(m,p)\S a ,m) a \S b ,p-m) b , 

(40) 

where p = if S a — Sb is an integer, while p = ±1/2 if 
S a — Sb is a half integer. 

For S a = Sb — S and J± ^> J z , we have calculated the 
entanglement entropy of \Gc), 

£{\Gc)) - -^[Z(p)/(m,p)] 2 log 2S+1 [Z(p)/(m,p)] 2 . 

(41) 

It is evaluated that when S — » 00, £(\Gc)) ~ 1/2, which 
is very large. 

When J± and J z approach the boundary J z = J± > 
from the regime of \Gc), it approaches I^q — Sb,p), where 
p = if 5 Q — 5b is an integer while p = ±1/2 if S a — Sb is 
a half integer. When Jj_ and J z approach the boundary 
J z = —J± < from the regime of \Gc), it approaches 
e^\S a + S bl p). 



D. J ± <-\J Z \ 

The energy spectrum for J± < can be obtained by 
using U{J±,J Z ) = UH(-J±,J Z )U*, where C/ = e" 5 "". 
Therefore, in the regime Ji < — \J Z \, the energy spec- 
trum is also given by Eq. 

The ground state is thus 



\G D ) = U\G C ) 

= z (p)/^2\f( r n,p)\\Sa,m)\S b ,p- 



(42) 
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with max(— S a , S z — Sb) < m < min(S a , S z + Sb)- 

Obviously the entanglement entropy of \Gjj) is the 
same as that of \Gc), with J± reversing its sign. 

When Jl and J z approach the boundary J z = —J± > 
from the regime of |Gx>), it approaches <"|£ — 
Sb,p). When J± and J z approach the boundary J z = 
J± < from the regime of \Gd), it approaches \S a + 
S b ,p). 



E. The ground states on the four parameter 
boundaries 
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When J z = Jj_ > 0, the Hamiltonian (JTJ is "H = 
J z S a -S b , the degenerate ground states are \S a — Sb, S z ) = 
J2 m 9(S a - S b ,S z ,m)\S a ,m)\S b ,S z - m), with S z = 
S b —S a , ■ ■ ■ , S a —S b . Here g(S a —S b , S z , m) is the Clebsch- 
Gordan coefficient. 

When J z = —J± > 0, the degenerate ground 
states are e™ S «*\S a - S b ,S z ) = £ m (-l) m <?(S a - 
S b , S x ,m)\S a ,m)\Sb, S z -m), with S z = S b -S a , ■ ■■ ,S a - 
Sb- 

When J z = J± < 0, the ground states are \S a + 
S b , S z ) = Y. m 9(Sa + S b , S z ,m)\S a ,m)\S b ,S z - m), with 
S z = —S a — S b , ■ ■ ■ , S a + Sb- 

When J z = — J± < 0, the ground states are e' l7rSaz \S a + 
S b , S z ) = 52 m {—l) m g(S a + S b , S z ,m)\S a ,m)\S b , S z -m), 
with S z = -S a — Sb, ■ ■ ■ , S a + Sb- 

The boundaries are where quantum phase transition 
take place. We have known that the ground states \Ga), 
\Gb), \Gc), \Gd), in the four regimes discussed in pre- 
vious subsections, depend on the values of J z and J±. 
Starting as a ground state in one of these regimes (see 
FIG. [8]), when J z and J± adiabatically approach each 
boundary regime, the ground state always approaches 
one of the degenerate ground states on the boundary. 
In entering the other regime across the boundary, the 
ground state restarts from another one of the degenerate 
ground states on the boundary. 



FIG. 9: (Color online) The entanglement entropy of the 
ground state as a function of I/77 == J z /Jj_. The filled 
symbols describe the analytical results. The empty symbols 
describe the numerical results. 



glected is regime B, i.e. 1/rj < — 1, as the ground state 
is exactly \S a , S a )\S b , S b ) or \S a , -S a )\S b , -S b ), without 
entanglement. Figure |H] clearly indicates excellent fitting 
between the analytical results in this section and the nu- 
merical results. 

Excellent fitting between our analytical results and the 
numerical results are also obtained for excited states. We 
have calculated the lowest energy states of different val- 
ues of S z for S a = 12000 and S b = 10000. Figure [TU] 
shows the the regime I/77 > — 1, i.e., regimes C and A, 
while FigfTTI shows the regime I/77 < —1, i.e., regime B. 
The reason for this separation is that the low-energy ex- 
cited state in regime B is with large magnitudes of S z , 
while those in regimes C and A are with small magni- 
tudes of S z . 

In conclusion, our analytical results fit the numerical 
results very well. 



IX. SUMMARY 



F. Comparison with the numerical results 

As each eigenstatc for J± < can be obtained by act- 
ing e mSaz on an eigenstate for J± = \J±\, we only need 
to consider the half of the parameter space with J± > 0. 

In this half parameter space, We have calculated the 
dependence of the entanglement on I/77 = J z /J±, using 
the ground states analytically obtained above in regimes 
A, C, and B. We compare these analytical results with 
the numerical results. The reason of choosing l/?y rather 
than 77 is because J z = in the middle of the half pa- 
rameter space. In this half parameter space, regime B is 
1/?/ < —1, regime C is — 1 < I/77 < — 1, while regime A 
is 1/n > 1. 

Figure [5] shows the the entanglement in the ground 
states for different values of S a = S b for I/77 > — 1. Ne- 



In this paper, we considered a binary mixture of two 
species of pseudo-spin- i atoms with interspecies spin ex- 
change in the absence of an external potential, and ex- 
tended the study of its ground states to the whole param- 
eter space of the two effective spin coupling strengths. 
Meanwhile, this provides a model of studying the rela- 
tion between the classical model and quantum ground 
states. 

We first analyzed the corresponding classical Hamilto- 
nian. We found the fixed points of the classical dynamics, 
and discussed their stability situation both analytically 
and numerically. The bifurcations were discussed. 

The classical evolution can be reproduced in quantum 
dynamics if starting from an initial state which is dis- 
entangled between the two species, as we have demon- 
strated. 

In the case that the atom numbers of the two species 
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FIG. 10: (Color online) The entanglement entropy of the 
lowest energy excited states for different values of S z under 
S a = 12000 and S b = 10000, as a function l/jj = J z /J±. Here 
we show the regime 77 > —1, in which the low-energy excited 
states are with small magnitudes of Sg. The filled symbols 
describe the analytical results. The empty symbols describe 
the numerical results. 
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FIG. 11: (Color online) The entanglement entropy of the 
quantum states as a function 1/rj = Jj_/J 2 , for different val- 
ues of S z under S a = 12000 and 5*6 = 10000 in regime B 
(I/?? < —1), where the low energy excited states have large 
magnitudes of S z . The filled symbols describe the analytical 
results. The empty symbols describe the numerical results. 



are equal, we confirmed in our system the previous claim 
that a classical fixed point bifurcation corresponds to 
maximal entanglement in the quantum ground state. 
Moreover, we find the result that when the two atom 
numbers are unequal, the entanglement of the quantum 
ground state at the parameter point of the bifurcation is 
not maximal, while the state corresponding to the fixed 
point that bifurcates indeed possesses maximal entangle- 
ment at that parameter point. 

A quantum ground state can be regarded as the clas- 
sical ground state with quantum fluctuations. This per- 
spective leads to solutions of the ground states in all pa- 
rameter regimes, by obtaining an effective Hamiltonian 
near each classical ground state. Using entanglement en- 



tropy as the quantity characterizing the ground states, we 
find that the analytical results fit the numerical results 
very well. We have made many detailed discussions. 

Our work establishes EBEC as a system manifesting 
connections between classical dynamics and quantum be- 
havior. 
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Appendix: Classical fixed points 

For our problem, the desired Lyapunov function will 
always be found by defining 



C = 71% + j2Jz(Saz + S bz 



(A.I) 



It is 



where 71 and 72 are suitably chosen coefficients, 
clear that jfc = 0. 

We find the following fixed points speci- 
fied by the values of S a = S Q n a , where 
n Q = (sin a cos <Pa , sin a sin (p a , cos a ) , < 6 a < 
ft, < (f a < 2n, a = a, b. 

(1) n a = — rib = (0,0, ±1); that is, one spin is 
parallel to the z direction, the other is antiparallel to 
the z direction. At these two point the eigenvalues 

of J arc fix = 



/i 2 = 0, H34 = ± 



C3-C4 



M5,6 



'^±Ci, where Ca = 2J 2 ± S a S b - J 2 z {S 2 a + S 2 ), C 4 = 

Jz(S a - S b )y/J*(S a + S h ) 2 - UJS^. If wc define m = 

5 ( \/W ' w ^ en 1 ^ z I ^ [ Jj_ 1 3 there are eigenval- 

ues with positive real part, and these two fixed points are 
unstable. Otherwise, the stabilization cannot be judged 
by the eigenvalues. One finds the Lyapunov function 
C = % — lf(S az + Sbz) 2 , which is minimal for J z > (or 
maximal for J z < 0) at each of these two fixed points in 
the parameter region 771 1 J z \ > | Jj_ | , where each of these 
two fixed points is thus stable. 

(2) n Q = —rib = (cos^s, sincp, 0), where < <p < 2n. 
The two spins are antiparallel and both are on the x — y 
plane. At this point, the eigenvalues of J are [i± = 1x2 = 
^ = /J4 = 0, M5, 6 = ±^/2J ± J z S a S b - J'i(S 2 + S 2 ). If 
we define 772 = 1^1% , when < J± < rj 2 Jz or rj 2 J z < 



J± < 0, some eigenvalues have positive real part, hence 
this fixed point is unstable. When J± > 772 J 2 > or 
J± < r]2J z < 0, one finds the Lyapunov function C = 
TL + y2Jz(S az + Sf, z ) 2 , which is minimal at the fixed point 
as 72 —J- 00. When J Z J± < 0, one finds that C = —H + 
l2Jz{S az -\-Sb z ) 2 is minimal at the fixed point as 72 — > 00. 
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TABLE I: Stable regimes of the fixed points. 771 = \ + \J§^ 



1 ~ s^+si ■ 



No. Fixed Points 

1 n a = -n b = (0, 0, ±1) 

2 n a — —Jib — (cos tp, sin ip, 0) 

3 n a = n b = (cos 93, sin <p, 0) 

4 n a = n 6 = (0,0, ±1) 

5 n a = -rif, 

6 n a = n 6 

7 n a = (sin cos ip, sin sin tp, cos 0) , 
nt = (sin 9 cos </?, sin 9 sin 99, — cos 6) 

8 n a = (sin 9 cos 75, sin 9 sin <p, cos 0) , 

n 6 = (— sin 6> cos 95, — sin # sin ip, cos 5) 



Stable regions 

m\M > \J±\ 

J± > 0,Jj_ > 772 Jz Or Jx < 0, J± < 772 Jz 
J±_ > 0, Jx > —J)lJz Or Ji < 0, Jx < — 772 

All 

Jx = Jz 

JX = Jz 

Jx = -J, 

Jx = — Jz 



Thus this fixed point is stable if Jj_ > 0,Jj_ > r\iJz or 
Jx < 0, Jx < ?72J 2 - 

(3) n a = rib = (cos sin (y9, 0); that is, the two spins 
are parallel and on the x — y plane. At this point, the 
eigenvalues of J are \i\ = [±i = [A3 = [14 = 0, /Lts^ = 
±y/-fl(S* + S£) - 2J A _J z S a S b . When - m J z < J ± < 
or < Jj_ < — f]iJz, some eigenvalues have positive 
real parts, hence the fixed point is unstable. For J±J Z > 
0, one finds £ = —TL + 72 J z {S a z + Sfcz) 2 is minimal at 
the fixed point as 72 — > 00. For J±J Z < one finds 
£ = H + 72 J z {S az + Sbz) 2 , which is minimal at the fixed 
point as 72 — > 00. Therefore the fixed point is stable 
when Jj_ > 0, Jx > — J72 Jz or J± < 0, J± < —r\iJz- 

(4) n = lib = (0,0, ±1); that is, the two spins 
are both parallel or antiparallel to the z direction. At 
each of these two fixed points, the eigenvalues of J 

are ^ = fx 2 = 0, /i 3 , 4 = ±y^^,/j 5 , 6 = ±y^±^; 
here Ci = -2JjS a S b - J 2 z {S 2 a + S 2 ), ( 2 = J z {S a + 
Sb)-\/J%(S a - Sb) 2 + AJ'^SaSb. The stabilization cannot 
be judged by the eigenvalues. But one finds the Lya- 
punov function £ = —(S az + Sbz) 2 , which is minimal at 
the fixed point. Hence these two fixed points are always 
stable. 

(5) In case J± = J z , the solution n a = n b with any 
possible is a fixed point; that is, the two spins are always 



parallel. The Lyapunov function £ = — S a ■ S b is minimal 
here, thus this fixed point is stable. 

(6) In case J± = J z , the solution n a = — rib with any 
possible is a fixed point; that is, the two spins are al- 
ways antiparallel. The Lyapunov function £ = S a ■ Sb is 
minimal here, thus this fixed point is stable. 

(7) In case J± = —J z , n a = 
(sin 8 cos ip, sin# sin tp, cos 6) while rif, = 
(sin 8 cos ip, sin# sin ip, — cos 9) is a fixed point; that 
is, the z components of the two spins arc opposite. 
One finds a Lyapunov function £ = — S a ■ S' b , where 
S b = (S'b sin 6 cosy, Sb sin 6 sin ip, Sb cos 6), which is 
minimal at this fixed point. Thus this fixed point is 
stable. 

(8) In case J± = — J z , n a = 
(sin 8 cos ip, sin# sin ip, cos 9) while rif, = 
(— sin 8 cos tp, — sin 8 sin tp, cos 8) is a fixed point; 
that is, the x and y components of the two spins are 
opposite. One finds a Lyapunov function £ = S Q • S fc ', 
where S' b ' = {—Sb sin 8 cos ip, —Sb sin 8 sin ip, —S b cos 8), 
which is minimal at this fixed point. Thus this fixed 
point is stable. 

All of the fixed points and their stable regimes are 
listed in Table |U 
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